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Abstract: We study boundary scattering in the </> 4 model on a half-line with a one- 
parameter family of Neumann-type boundary conditions. A rich variety of phenomena is 
observed, which extends previously-studied behaviour on the full line to include regimes of 
near-elastic scattering, the restoration of a missing scattering window, and the creation of 
a kink or oscillon through the collision-induced decay of a metastable boundary state. We 
also study the decay of the vibrational boundary mode, and explore different scenarios for 
its relaxation and for the creation of kinks. 
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1 Introduction 

Systems with boundaries, defects and impurities have been intensively studied in statistical 
physics and field theory, both at the classical and the quantum levels. Often the key 
physics of the model can be captured, possibly after dimensional reduction, by a simple 1+1 
dimensional field theory on a half line. Examples include the Kondo problem [1], fluxon 
propagation in long Josephson junctions [2], the XXZ model with boundary magnetic field 
[3], an impurity in an interacting electron gas [4], the sine-Gordon [5] and Toda [6] models, 
monopole catalysis [7], the Luttinger liquid [8], and a toy model motivated by M-theory [9]. 

Especially since the work of Ghoshal and Zamolodchikov [5], there has been great 
interest in boundary conditions compatible with bulk integrability, and many such models 
turn out to be of direct physical interest. However less attention has been paid to the 
equally if not more physically-relevant cases of non-integrable boundary systems, even at 
the classical level. This is perhaps a shame, as it is now known that non-integrable classical 
field theories, even in 1+1 dimensions, can exhibit remarkably rich patterns of behaviour 
not seen in their integrable counterparts [10-14]. 
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In this paper we examine the 0 4 theory in 1+1 dimensions, restricted to a half line 
by a simple Neumann-type ‘magnetic field’ boundary condition. (The sine-Gordon model 
with a non-integrable boundary was recently investigated in [15].) The 0 4 theory on a 
full line is similar to the sine-Gordon model in that both support topological kinks and 
antikinks; the 0 4 theory also has an intriguing and still not fully-understood counterpart of 
the sine-Gordon breather, the oscillon [16]. We chose the magnetic field boundary condition 
in part because of its simplicity, and in part because the scattering of kinks against such 
a boundary provides a natural deformation of the full-line scattering problems which are 
already known to exhibit intricate patterns of resonant scattering [10-13]. In some regimes 
our results do indeed resemble the pattern of scattering windows observed in kink-antikink 
collisions on the full line, while in others we find novel phenomena including a new type of 
‘sharp-edged’ scattering window. Even though the theory is not integrable, it turns out to 
be possible to give an accurate analytical description of some aspects of this behaviour. We 
complement these studies with an investigation of the decay of the vibrational boundary 
mode through nonlinear couplings to scattering states, and of the creation of kinks by 
an excited boundary. An interesting feature of the boundary mode decay, discussed in 
section 8 , is that with suitable initial conditions a period of relatively slow decay can be 
followed by a sudden burst of radiation from the boundary as a new decay channel opens. 

While this paper is self-contained, we have also made a number of short movies to illus¬ 
trate aspects of the discussion, which can be found at [17]. After a brief explanation of the 
numerical methods used to obtain our plots in appendix A, these are listed in appendix B. 

2 The model 

We consider a rescaled 0 4 theory with vacua c/) v G {—1, +1} on the left half-line —oo<x< 
0. The bulk energy and Lagrangian densities are £ = T + V and C — T — V respectively, 
where 

T = + 2 and V = \<t? x + + 2 - l ) 2 . ( 2 . 1 ) 

The static full-line kink and antikink, (/>k{x) = tanh(x — xo) and (/>x( x ) = ~ ( I ) k( x )^ have 
rest mass M — 4/3 and interpolate between the two vacua. Including a boundary energy 
—H<j) o, where 0o = 0(0, t) and H can be interpreted as a boundary magnetic field, yields 
the Neumann-type boundary condition 0^(0, t) — H at x = 0. 

For 0 < H < 1 there are four static solutions to the equations of motion, shown in 
figure 1. Two of them, 0 \{x) — tanh(x — Xo) and 02 (x) — tanh(x + Xo) with Xo = 
cosh+l /yfiH\), are restrictions of regular full-line kinks to the half-line, while the other 
two, 0 s{x) — — coth(x — Xi) and 04 (x) = — coth(x + Xi) with X\ — sinh _ 1 (l/+ / |i7|) 
are irregular on the full line. On the half line, 0 3 is non-singular and corresponds to the 
absolute minimum of the energy, while 01 is metastable, and 02 is the unstable saddle-point 
between 0 3 and 0Their energies can be found by rewriting £[0] = J_ Vcte — £0 0 in 
Bogomolnyi form as 

0 

m = 5 /+ ± (<+!))+ T [+++ - Hcf> 0 . ( 2 . 2 ) 
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Figure 1. Static solutions for H = 1/2. 

Since 0i and 02 satisfy <fi x = 1 — 0 2 we have 0i(O) = —y/l—H, 02 ( 0 ) = y/l—H ; while 
(0 s) x = 03 — 1 and so 03(0) = y/l+H. Taking the upper and lower signs in (2.2) as 
appropriate, 

e[M = I - la-#) 3/2 , em = | + fa-#) 3/2 , 

^ 3 ] = |-i(l+^) 3/2 . (2.3) 

As H increases through 1 , 0 1 merges with 02 and disappears, leaving 0 3 as the only 
static solution for H > 1. For H < 0 the story is the same, with 0 and H negated 
throughout, so the physically-relevant solutions are (j)i(x) := — 0z(x), i — 1... 3. 

3 Numerical results 

We took initial conditions corresponding to an antikink at xo = —10 travelling towards 
the boundary with velocity > 0. (We found the setup with an incident antikink easier 
to visualise, but our results apply equally to kink-boundary collisions on negating 0 and 
H .) Thus the initial profile was 0(x,O) = 0i(x) — tanh(y(x — xo)) + 1 for H >0 and 
0(x, 0) = 03(x) — tanh( 7 (x — xo)) + 1 for H < 0, where 7 = 1/yjl — vf. Our real interest 
was in the problem with the initial antikink infinitely far from the boundary, but the rapid 
decay of the antikink-boundary force (4.1), calculated below, meant that error in taking 
xq — —10 was small. 

To solve the system numerically, we restricted it to an interval of length L, with the 
Neumann boundary condition imposed at x = 0 and a Dirichlet condition at x = —L. 
(Since we took run times such that radiation did not have time to reflect from the extra 
boundary and return, the boundary condition at x — —L was anyway irrelevant.) We used 
a 4 th order finite-difference method, explained in more detail in appendix A, on a grid of 
N = 1024 nodes with L = 100, so the spatial step was 5x ~ 0.1, and a 6 th -order symplectic 
integrator for the time stepping function, with time step St = 0.04. Selected runs were 
repeated with other values of xq, T, N and St to check the stability of our results. 
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Figure 2. Final antikink velocities as functions of initial velocities. The dashed line indicates the 
result for a purely elastic collision. In the fifth plot, a kink can also be produced: its velocity is 
shown in red. The horizontal dotted lines in plots a and b show the relevant values of v cr (H), as 
given by equation (4.2) below. 


Our simulations revealed a rich picture, aspects of which are summarised in figures 2 
and 3. For all (iF, vi) pairs with H < H c « 0.6, the antikink either reflects off the boundary 
with some velocity uj, or becomes stuck to it - corresponding to Vf = 0 - to form a ‘bound¬ 
ary oscillon’. This latter configuration oscillates with a large (of order one) amplitude, and 
a below-bulk-threshold basic frequency. Just like the bulk oscillon (which it becomes in the 
limit H 0), it then decays very slowly into radiation. At H = 0 (figure 2c) the plot of 
\vf\ as a function of Vi reproduces the well-known structure of resonant scattering windows 
in KK collisions on a full line [10-12]. For negative values of H (figures 2a and b) new 
features emerge. For small, the antikink is reflected elastically from the boundary with 
very little radiation. As increases above an iF-dependent critical value v cr: the antikink 
is trapped by the boundary, leaving only radiation in the final state. Increasing V{ further, 
scattering windows begin to open, until V{ exceeds an upper critical value and the antikink 
again always escapes. If the antikink does escape, its speed \vf\ is always larger than some 
minimal value very slightly lower than v cr: so (in contrast to the full-line situation) vj is 
a discontinuous function of giving the windows the sharp edges mentioned in the intro- 
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Figure 3. A ‘phase diagram’ of antikink-boundary collisions. The plot shows the value of the field 
at x = 0 a time tf = \xo\/vi + 100 after the start of the simulation, as function of the boundary 
magnetic field H and the initial velocity V{. 


duction. For small positive values of H (figure 2d), Vf is instead a continuous function of 
Vi, the sequence of windows for H — 0 shifting towards lower velocities while preserving its 
general structure. Finally, for H > H c (figure 2e) other new phenomena arise which have 
no counterparts in the full-line theory; these will be discussed further in later sections. 

4 Kink-boundary forces and the location of the low-velocity window 

To understand the novel window of near-elastic scattering at low initial velocities when H is 
negative, seen in figures 2a and 2b, we start by evaluating the static force between a single 
antikink and the boundary. Placing the antikink at x = xo < 0, we add a possibly-singular 
‘image’ kink at x\ > 0 in such a way that the combined configuration satisfies the boundary 
condition at x — 0. From the standard full-line result, the force on the antikink from the 
image kink is equal to 32e~ 2 ( Xl ~ Xo \ or minus this if the image kink is singular. For \H\ <C 1 
and \xo\ 1 we find that the boundary condition requires e~ 2xi — \H + e 2x °, so 

F = 32 [\H + e 2xo ) e 2xo . (4.1) 
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Figure 4. Example collisions for H = —0.4 (left column) and H = 0.9 (right column), illustrating 
various scattering scenarios. At negative H\ (a) elastic recoil for low impact velocity; (b) saddle 
point production at the critical velocity, the antikink finishing on the top of the potential barrier; 
(c) single bounce with the antikink escaping back over the barrier. At positive H : (d) single 
bounce with the excitation of the H > 0 boundary mode; (e) kink production via collision-induced 
boundary decay; (f) bulk oscillon production. See also movies M01 - M07 and Mil of appendix B. 


For H < 0 the force is repulsive far from the boundary, only becoming attractive nearer in. 
When xo — ^log(— \H) 1 x\ — oo and the force vanishes, the antikink-kink configuration 
reducing to the unstable static solution fa. 

Now consider, again for H < 0 , an antikink moving towards the boundary. If its 
velocity is small, then it won’t have sufficient energy to overcome the initially-repulsive 
force, and it will be reflected without ever coming close to x — 0 , and without significantly 
exciting any other modes; this behaviour is illustrated in figure 4a. Increasing at some 
critical value v cr the energy will be just enough reach the top of the potential barrier and 
create the static saddle-point configuration fa, as shown in figure 4b. The value of v cr can 
be deduced on energetic grounds: the initial energy is |(1 — v 2 r )~ 1 / 2 + E[fa], while the final 
energy is E[fa\ = § + |(1 +H) 3 / 2 . Equating the two, 

v cr (H) = sj 1 - 4((l+tf) 3/2 + (1-#) 3/ T 2 • (4.2) 

If Vi is just larger than u cr , the antikink can overcome the potential barrier and approach the 
boundary; energy is then lost to other modes and so it is unable to return, and is trapped at 
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the boundary. Thus v cr (H ) marks the upper limit of the windows of almost-perfectly-elastic 
scattering seen in figures 2a and 2b, and the lower edge of the ‘fractal tongue’ occupying the 
left half of figure 3. The curve Vi = v cr (H ) is included in figure 3; it matches our numerical 
results remarkably well. Indeed, it can be seen from figure 7 below that the maximum 
error is of the order of 0.5%, which is rather small given that radiation was ignored in the 
derivation. Similar arguments show that, within this approximation, v cr is the smallest 
possible speed for any escaping antikink, explaining the sharp (discontinuous) edges of all 
windows when H < 0. 

5 The boundary mode 

Next we consider the perturbative sector of the model, that is solutions of the form 0(x, t) = 
(/) s (x)+r](x, t) where (j) s {x) is a static solution to the equations of motion and rj(x , t) is small. 
The full-line theory has a continuum of small linear perturbations about each vacuum with 
mass mn — 2, while a static kink </)k(x) — tanh(x — Xq) has two discrete normalizable 
modes - the translational mode, and a vibrational mode with frequency uo\ — >/3 - and a 
continuum of above-threshold states rj(x,t) = e lujt r]k(x) where uo 2 — 4 + k 2 and [19] 

rjk(x) = e ~ l K x ~ x o) (_j — A: 2 + 3ifctanh(x — Xq) + 3tanh 2 (x — Xq)) . (5.1) 

Turning now to the half-line theory, we can regard (/)k(x) instead as the static half-line 
solution 0i (x) to the boundary theory with 0 < H < 1 and 0(—oo) = — 1. Its linear 
perturbations must now satisfy d x r\[x ) = 0 at x — 0. Setting k = ift this yields 

ft 3 — 30oft 2 + (60o — 4) ft — 60o + 60o = 0 (5.2) 

where 0o = 0i(O) = —y/1 — H and the frequency ujb of the corresponding boundary mode 
satisfies uj 2 b = 4 — ft 2 . The solutions of (5.2) for both negative and positive values of 0o are 
shown on the left-hand plot of figure 5; note that only solutions with ft > 0 can give rise to 
localised modes, and of these, ft must be less than 2 for ujb to be real and the mode stable. 
We will denote the corresponding normalised profile function as t]b{x) := rji^x)/ : 
where rji K (x) is given by (5.1) with k = m. 

For 0 < H < 1, we have — 1 < 0o < 0 and (5.2) has just one positive solution ft, which 
satisfies ft < 2: this is the single vibrational mode, localised near to the boundary. The 
linear perturbations of 02 (x), the saddle-point solution, are also described by (5.2), but 
now with 0o = 02 ( 0 ) = +y/l H. For these cases (5.2) has two positive solutions but one 
is larger than 2: this is the unstable mode of 02 (x). Finally, for H < 0, the continuation 
of (5.2) to 0o < —1 governs the spectrum of fluctuations about 03(x), the H < 0 vacuum 
in the 0(—oo) = —1 sector. There are no positive solutions in this regime and hence no 
vibrational modes of the boundary for H <0. The right-hand plot of figure 5 summarises 
the situation, plotting the images of the positive-ft parts of the curves shown on the left 
under the mapping (0o, ft) —> = ( 1 —0q, 4—ft 2 ). The grey dashed parts of the curve 

visible for H < 0 are included for completeness but do not describe vibrational modes of 
physical solutions - they correspond to ‘perturbations’ of the singular solution 04 (x). 
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Figure 5. Linearised boundary mode analysis: on the left, the solutions of equation (5.2) as a 
function of 0o = 0(0); on the right, the frequencies of localised boundary modes as a function of H. 
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Figure 6. Power spectra at the boundary after a collision with vi = 0.5, for H = —0.1 (upper) 
and H = 0.3 (lower). 

These findings are confirmed by our numerical results. Figure 6 shows the Fourier 
transforms of 0(0, t) for 30 < t < 3030, for antikink-boundary collisions with initial velocity 
Vi = 0.5, and H = —0.1 and 0.3. The final velocity Vf of the reflected antikink is —0.382596 
for H = —0.1 and —0.454014 for H = 0.3, so in both cases translational energy is lost to 
other modes during the collision. 

For H — —0.1, the boundary does not have an internal mode, and only radiative 
modes with frequencies near to 2, the mass threshold, remain near to the boundary. The 
internal mode of the reflected antikink has frequency c ji, but this mode cannot be observed 
at the boundary since it is exponentially suppressed there. However nonlinear couplings 






































with other excitations create waves with frequencies at above-threshold multiples of uj\ 
[20], which can propagate back to the boundary. The upper plot of figure 6 shows peaks at 
Q.i — 2 and D 2 = D(2cji), where 0(cj) = 7 (cj + k(uj)vf) is the Doppler-shifted frequency of 
radiation emitted from the moving kink measured on the boundary. Higher harmonics at 
D 3 = 0(3cji) and D 4 = D(4c<;i) are also visible, along with combinations of the internal mode 
of the antikink and the lowest continuum mode such as D 5 = 0(2+uji) and Dg = 2+D(2 cji). 

Many of these modes are also present in the H — 0.3 spectrum shown in the lower plot of 
figure 6 , albeit at shifted locations because of the different final antikink velocity. However 
the plot is dominated by the internal boundary mode with frequency Dio — w b — 1.888459. 
The higher harmonics Du = 2c ob and D 12 — 3ub are also visible, while interactions between 
radiation from the outgoing antikink and the boundary mode lead to peaks at D 13 = 
ujb + D(2cji) and D 14 = uob ~ D(2ccq). 

6 The resonance mechanism in boundary scattering 

For small nonzero values of \H\, the resonant energy exchange mechanism governing scatter¬ 
ing in the bulk 0 4 model is changed in two ways in the boundary theory: (i) the attractive 
force acting on the antikink near to the boundary is modified, in particular becoming re¬ 
pulsive at greater distances when H is negative; (ii) after the initial impact, energy can 
be stored not only in the internal mode of the antikink, but also, for positive values of 
i7, in the boundary mode. These factors change the resonance condition for energy to 
be returned to the translational mode of the antikink on a subsequent impact after some 
integer number of oscillations of the antikink’s internal mode, shifting (and, for negative 
H , sharpening) the windows seen in figures 2 a-d. This return can happen after two, three 
or more bounces from the boundary, leading to a hierarchy of multibounce windows as in 
the full-line situation. Our numerical results suggest that for small positive values of H the 
contribution of the boundary mode in the resonant energy transfer is not significant. 

For larger values of \H\ other new features appear. For H < 0 the first is the resurrec¬ 
tion of a two-bounce window that was observed to be missing from the full-line scattering 
process by Campbell et al in [10]. Figure 2a includes a scattering window centred at 
Vi ~ 0.245 which is not the continuation of any of the windows seen in figures 2b-d; the 
same window can be seen in figure 3 running from (H,Vi) = (—0.135,0.202) to (H,Vi) = 
(—0.489, 0.417), and the top plot of figure 7, running from (i7, Vi~v cr ) = (—0.135, 0.084) to 
(i7, Vi) = (—0.489,0.023). The emergence of this window as H decreases below H ~ —0.136 
is shown in more detail in the left-hand set of plots of figure 8 . The nature of the new win¬ 
dow is made clear by the plots in figure 9, which shows 0(0, t) for V{ inside the first three 
two-bounce windows for H = —0.2, and also for H — 0, which is equivalent to the full-line 
case. The ‘wobbles’ between the large dips in such plots count the oscillations of the internal 
antikink mode between bounces [ 10 ]. As can be seen from the figure, the minimum num¬ 
ber of oscillations supporting antikink escape is one smaller for H — —0.2 than it was for 
H — 0, giving rise to the extra window. A complementary process of ‘window destruction’ 
occuring for H > 0 can be seen on a close examination of figure 3, and on the right-hand 
set of plots of figure 8 . Decreasing H further, we also observed interesting structures at the 
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Figure 7. Zoomed-in views of the region near the tip of the fractal tongue of figure 3. Note that 
the vertical axes show multiples of Vi — v cri where v cr = v cr (H ) is the theoretical upper limit of the 
near-elastic scattering window given by the formula (4.2). 
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tip of the fractal tongue, near to H = — 1, with resonance windows merging to give rise to 
a pattern of half-rings on the phase diagram. These are shown in the lower two plots in 
figure 7. 

For H > 0, the scattering can induce the metastable (j) i boundary to decay to 3, the 
true ground state, with the creation of an extra kink. This process, which has no analogue 
in the full-line theory, is visible in the extra red ‘kink’ line in figure 2e. The principal 
region of boundary decay occupies the solid red area on the right edge of figure 3, and is 
examined in more detail in figure 10. If the boundary mode is sufficiently strongly excited 
by the initial antikink impact, it behaves as an intermediate state prior to the escape of 
a kink from the boundary, analogous to the intermediate oscillon state in the process of 
KK pair production on the full line [18, 20, 21]. Depending on their relative velocities, 
the reflected antikink and the subsequently-emitted kink may appear separately in the final 
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Figure 8. Scans of 0(0, t ) for various values of H , showing the emergence (on the left) and 
destruction (on the right) of a two-bounce window as H moves away from zero. 



Figure 9. Plots of 0(0, t ) for Vi inside the first three two-bounce windows for H = —0.2 (left) and 
for H = 0 (right). 
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Figure 10. Antikink-boundary scattering at large H. (a): a zoomed-in view of figure 3 showing 
the value of the field at x = 0 a time tf = \xo\/vi + 100 after the start of the simulation; (b): the 
measured final velocity of the reflected antikink; (c): the measured final velocity of the emitted 
kink, if present, white being plotted otherwise; (d): the difference between these two velocities. 


state, or recombine to form a bulk oscillon. Such collisions lead themselves to a fractal-like 
structure with windows where the antikink and kink separate interspersed with regions of 
oscillon production, just as in the full-line theory (though with added complications due 
to interference with radiation from the boundary). Some of this structure can be seen 
in figure lOd, where the blue regions inside the zone of boundary decay show windows 
of antikink and kink separation, while the yellow regions correspond to the production 
of a bulk oscillon, and also in the movies Mil and M12. Spacetime plots of some of the 
relevant processes, for H — 0.90, are shown in the right panels of figure 4: scattering of 
the antikink with excitation of the boundary mode, but no kink production (d); production 
of a separated KK pair, with the boundary decaying to the true ground state (e); and 
recombination of the KK pair to form a bulk oscillon (f). A further intriguing feature 
of the region of boundary decay, clearly visible in figures 3 and 10, is the cusp-like nick, 
terminating at (i7, Vi) ~ (1, 0.365), which splits it into two disconnected parts. This appears 
to be associated with a velocity-dependent vanishing of the effective coupling between the 
incident antikink and the boundary mode. It would be very interesting to find an analytical 
understanding of this phenomenon, but we will leave this for future work. 
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7 Radiative decay of the boundary mode 

A significant feature of the </> 4 model is that its spectrum of perturbative oscillations around 
the static kink or antikink solutions contains an internal vibrational mode. If the amplitude 
of the excitation is small enough and nonlinear corrections can be neglected, this mode os¬ 
cillates with almost-constant amplitude A and frequency — v^3. For larger amplitudes 
nonlinearities start to play an important role. It has been shown [20] that the first anhar- 
monic correction to the internal mode oscillation results in the appearance of an outgoing 
wave with frequency 2^, which is above the mass threshold. The corresponding rate of 
radiative energy loss is dE/dt A 4 , causing the mode to decay. The resulting time depen¬ 
dence of the amplitude of the internal mode follows the law dA/dt ^ A 3 , where the explicit 
value of the proportionality constant can be found using a Green’s function technique [20]. 

For our boundary theory, we have observed a similar pattern in the decay of small- 
amplitude excitations of the boundary mode, but with a number of interesting new features. 
For small positive values of i7, the frequency ujb of the linearised boundary mode, as 
predicted by (5.2), satisfies 2 ujb > 2, and so the second harmonic of this mode is able to 
propagate in the bulk 1 . But as the boundary magnetic field H increases, the frequency of the 
boundary mode decreases, and when H > H 2 ~ 0.925, 2 ujb dips below 2 and the situation 
changes. The second harmonic can no longer propagate into the bulk, and this channel 
of radiative energy loss from the boundary is terminated. Only the next harmonic, which 
appears in the third order of the perturbation series, can be seen in the power spectrum. 
The radiation loss rate becomes dE/dt ~ A 6 and the decay rate is reduced to dA/dt A 5 . 

The situation changes again as H increases beyond H — Hs « 0.982, when 3 ujb 
falls below 2 and the third harmonic joins the second, trapped below the mass threshold. 
Theoretically, as H 1 and ujb 0 this pattern will repeat an infinite number of times, 
so that whenever ^j-j- < wb < f, the amplitude of the decaying mode should satisfy, to 
leading order, the equation dA/dt A 2n-y 

Figure 11 shows the behaviour of the field on the boundary and at x = —50, in the 
far field zone, with initial conditions (f)(x, 0) = + 0.05//^(x), 0) = 0 and H = 

0.90 < i? 2 - The power spectrum of the field on the boundary is dominated by boundary 
mode oscillating with the theoretically predicted frequency ojb = 1.08509. There are also 
two peaks at 2 ujb and 3 ujb- Since ujb < m, this lowest mode cannot propagate and indeed, 
there is no trace of it in the far field zone. The mode with the frequency 2 ujb is already 
in the scattering spectrum, so this mode does propagate, causing the energy loss from the 
boundary mode, as seen in figure 11a. 

The picture is different when H = 0.94 > H 2 (see figure 12). The mode with frequency 
lob — 0.93643 still dominates the power spectrum of the boundary excitations, but its decay 
is much slower, reflecting the fact that the mode 2 ujb is now below the mass threshold and 
cannot propagate into the bulk. As can be seen from the power spectrum in the far field 
zone plotted in figure 12d, the radiation is much less than in the previous case. There are 
two dominant frequencies, uo = 2 and cj = 3cjb- The presence of the peak at 3 wb is natural, 
since this is the first harmonic above the mass threshold. The peak at uo = 2 originates from 

1 Recall that m — 2 is the mass threshold for the bulk theory. 
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Figure 11. Evolution of the boundary mode for H = 0.90 < H 2 , with initial conditions 0(x, 0) = 
+ 0.0 5t]b(x), 0) = 0 : a) the amplitude of the boundary mode as a function of time; 

b) its power spectrum, found by taking the Fourier transform of 0(0, t) for 0 < t < 1600; c) the 
values of r)(x,t) := 0(x,£) — 0 \(x) at x = —50; d) its power spectrum, taken from rj(— 50, t) for 
200 < t < 1600. In plots a) and c), solid non-transparent lines join local extrema of the measured 
field at the given position. 




Figure 12. The same sequence of plots as in figure 11, but now for H = 0.94 > H 2 . 


near-threshold bulk modes, excited by the initial conditions via the nonlinearities, which 
disperse only slowly away from the boundary [22]. 
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Figure 13. Maximal radiation amplitude l^maxO^)!? where r)(x,t) = (/>(x,t) — (j)\{x) is the deviation 
of the field from its static value, at x = —50 for the kicked initial conditions </>(x,0) = </>i(x), 
0t(x,O) = AqujbVb{x) as a function of H , for three different values of A 0 . 



Impetus A 0 

Figure 14. A log-log plot of the maximal radiation amplitude at x = — 50 for the kicked initial 
conditions </>(x, 0) = 0i(x), </>t(x,0) = Aqujb^b{x) as a function of A$. Note that for small Aq and 
H < H 2 the power law decay rate is universal. 


Another test of the scenario is to consider the radiation from a “kicked” boundary 
initial condition </>(x, 0) = </>i(x), 4>t(x, 0) = A^ujbtib{x) in the far field zone. Our numerical 
results for this case are presented in figures 13 and 14. 

Figure 13 shows the independence of maximal amplitude of the field measured at 
x = — 50, far away from the boundary, for three small values of the initial impetus Aq given 
to the boundary mode. Note that the radiation amplitude drops sharply when H crosses 
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H _2 and i? 3 , as predicted by our general considerations. 

Figure 14 shows a log-log plot of the dependence of the maximal amplitude at x = —50 
on Aq. For small values of Aq and H < all curves have the same slope, fitting the 
expected ~ Aq dependence. The curve for H — 0.95 shows a significant reduction in the 
radiation amplitude, reflecting the loss of a decay channel as H passes i? 2 - However its slope 
for small values of Aq appears to be relatively unchanged from that of the previous curves, 
even though our previous considerations based on the propagation of the third harmonic 
would suggest an ~ Ajj dependence. It may be that slow (near-threshold) bulk modes, 
visible in figure 12d in the peak at uo — 2, are obscuring the effect we are looking for. It 
is possible that this could be tested by waiting significantly longer before measuring the 
radiation, to allow the slow modes to die away, but a more-detailed study would be needed 
to draw a clear conclusion. 

Another interesting feature visible on each curve is that as Ao reaches some (curve- 
dependent) critical value, the radiation flux suddenly dips. As will be discussed in the next 
section, this effect is associated with the nonlinear effect of the reduction in the frequency 
of the boundary mode with increasing amplitude. 

Finally, for even larger values of the intial impetus we can see a large increase of the 
amplitude of the field in the far zone. This is a signature of a non-perturbative effect, the 
excitation at the boundary becoming strong enough to destabilise it completely, with the 
emission of a kink into the bulk flipping the field there into the other vacuum. Some further 
observations concerning this phenomenon are reported in section 9 below. 

8 Higher-order nonlinear effects and amplitude-dependent decay rates 

In the last section we principally considered boundary mode decay in the small-amplitude 
regime where the boundary mode itself could be treated linearly. For larger amplitudes the 
frequency of the mode’s oscillation is lowered, just as in the case of an anharmonic oscillator 
or the simple pendulum. Numerical simulations of the oscillations of a full-line kink [20] 
also exhibit this behaviour, which is typical for many nonlinear systems. 

In the evaluation of the critical values H n above, we implicitly assumed that the ampli¬ 
tude of the excitation was small, so that its frequency was that predicted by the linearised 
equations. However for larger amplitudes, given the amplitude-dependent frequency re¬ 
duction just described, it is possible that even for H < i? 2 , the actual frequency of the 
boundary mode, Cob , will be lower than m/2. Then the decay rate will be slower than that 
observed for smaller amplitudes, since the second harmonic will not couple directly with 
any propagating bulk modes. However, the amplitude of the boundary mode will decrease 
with time due to the outgoing radiation, causing its frequency to grow. Provided H < i? 2 , 
once the amplitude has decreased far enough, the second harmonic will enter the scatter¬ 
ing spectrum. In such a case we can expect to observe an intriguing phenomenon: while 
initially the radiation flux from the boundary is relatively small and the decay rate rather 
slow, after some time there will be a sudden increase of the radiation flux and a switch to 
a much faster decay rate. 
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Figure 15. Evolution of the boundary mode for H = 0.8393 < H 2 and large initial amplitude 
A 0 = 0.3. Plots a) and c) show the time evolution of the amplitudes of the boundary mode and 
radiation field respectively, from the values of the field at x = 0 and x = — 50. Plots b) and d) show 
power spectra at the positions x = 0 and x = —50. The black lines show the power spectra from 
0(0, t) for times 0 < t < 1600 (plot b)) and 0(—50, t) for times 200 <t< 1600 (plot d)). The purple 
and green lines and filled areas show the power spectra for times before and after the transition 
(0 < t < 750 and 750 <t< 1600 for x = 0; 200 <t< 800 and 800 <t< 1600 for x = —50). 


Numerical work confirms that this effect really exists, as can be seen in figure 15 and 
movie M13, which show the decay of the amplitude of the boundary mode. For about the 
first 750 units of time the amplitude changes very slowly, albeit with a small modulation, 
after which there is a sudden transition to a much more rapid decay. In the far field zone 
this effect can be observed as a sudden jump of the radiation flux, by about one order of 
magnitude. 

In the power spectrum plotted in figure 15d one can clearly see a large peak just below 
uj — 1, which is the initial frequency of the mode. While the amplitude slowly decreases 
the frequency grows until it crosses 1, after which point the decay runs much faster. We 
can also see a drift of the frequency up to ujb — 1.24666. 

This slow-then-fast behaviour is reminiscent of higher-dimensional oscillon decay. In 
[23-25] it was observed that oscillons in two and three spatial dimensions lose their energy 
very slowly for tens of thousands of oscillations until they reach some critical frequency, 
above which they quickly decay to the vacuum. 
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9 Creation of kinks from an excited boundary 

The final phenomenon we investigated was the creation of kinks from the metastable bound¬ 
ary. We previously observed that this could be induced in certain scattering processes at 
large H. To view it in isolation, we instead excited the boundary mode directly, taking 
initial conditions of the two types (“stretched” and “kicked”) used earlier. First, we used 
initial condition 

0) = 4>i(x) + A 0 r] B (x), 0) = 0 (9.1) 

with t]b{x) the boundary profile for the linearised prolem, as an approximation to the 
boundary mode at its largest deviation from equilibrium; and second, we took 

4>{x,0) = (j>i(x), 0) = A 0 u B r]B(x) (9.2) 

representing the “kicked” boundary. As before, we normalized the profile of the boundary 
mode in such a way that tjb ( 0 ) = 1. 

In both cases, if Aq is taken to be sufficiently small, the boundary oscillates with 
frequency ujb and the amplitude Aq. However, as Ao becomes larger, the nonlinear processes 
discussed above start to play a significant role and further, as the initial energy of the excited 
mode becomes sufficient, outgoing kinks can be observed in the far zone, as seen in figures 16 
and 17. 

Note that for large H 1 the boundary mode profile resembles the difference between 
the unstable boundary solution (a saddle point of the energy) and the stable boundary 
solution: 

tanh(x + Xq) — tanh(x — Xq) 

2 tanh(Xo) 

Therefore the boundary mode, with appropriate amplitude, being added to the static 
boundary solution yields the unstable boundary. When the solution crosses the saddle point 
of energy it decays into another static solution with an additional kink is emitted from the 
boundary. 

Therefore the critical value of the amplitude of the boundary mode for the production 
of the kinks is: 

A C rit ~ 02 ~ 01 = 2Vl — H. (9.4) 

This critical amplitude is in very good agreement with the first type of the initial conditions 
for positive values of Aq. Only for very small values of Aq is there a symmetry Aq —>► — Ao. 
For larger Ao > 0 the excitation have less energy than the excitation for — Ao. Therefore 
the critical line for kink creation, from the left side of the plot, is much closer to the centre 
(A) = 0). 

For initial conditions of the second type, the energy for Ao and — Aq is exactly the same 
and therefore the plots look much more symmetric. For H — » 1 the critical amplitude is 
almost exactly A cr u = y/1 — H , half as big as in the first case. 


(9.3) 


Vb(x) « 
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Figure 16. The field value on the boundary at time tf = 50 for initial conditions 0) = 
c/)b(x) + Aqtib{x), 0) = 0. The blue colour represents the region without kink creation. The 
dashed line corresponds to the function H = 1 — A 2 /A. 
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Figure 17. The field value on the boundary at time tf = 50 for initial conditions 0) = (/>b(x), 
0) = AoujbVb{x). The blue colour represents the region without kink creation. The dashed 
line corresponds to the function H = 1 — A 2 . 
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10 Conclusions 


Our investigations of the boundary </> 4 theory have shown that it offers a considerably richer 
variety of resonance phenomena than the bulk theory, within a setting where analytical 
progress can be made. Key features include the modification of the force leading to the 
sharpening of window boundaries and the new critical velocity v cr: the resurrection of the 
first ‘missing’ scattering window, the observation of the boundary oscillon, and the collision- 
induced decay of the metastable boundary vacuum for H near to 1. Much of our work has 
been numerical and many issues remain for further study, the most pressing being the 
development of a reliable moduli space approximation incorporating the boundary degrees 
of freedom (see [9] for some earlier work on this issue). This model is sufficiently simple 
that it should offer the ideal playground for the development of better analytical techniques 
for the understanding of more general nonintegrable field theories. 
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A Numerical methods 


In this Appendix we describe some details of the numerical methods used in our simulation. 
In our numerical code we used the following discretization: 

u n = 4>(—nh), n = 0...N. (A.l) 

To calculate spatial derivatives we used a fourth-order central difference scheme 


D 2 u n = 


1 


12 K 2 


( ^n—2 T 10^77,-1 30 U n T \6u n -\-i 7i ri _|_2) 


(A.2) 


for all points far enough from the boundary, n < 2. This scheme can be derived using 
Lagrange polynomial approximation: 


n-\-m 


n-\-m 


u{x ) = y u ^( x )’ = n ^ 


i—n—m 


j=n—m 


X — jh 

ih — jh 


(A.3) 


However for the two points closest to the boundary we have to use a different basis 

3 


i(x) = H£ 0 (x) + Ui$i(x), 


(A.4) 


i =0 
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where 


h>o(x) 

x rr x-jh 
(■ ih) 2 ih — jh %> 

(A.5) 




£ 0 (x) 

d +1 e*)n x -T 

(A.6) 

£ 0 (x) 

’•ir/ 

(A.7) 


Note that for points within the interpolation intervals 

= 5ij, I(xi) = 0, Z'(0) = 0, £' (0) = 1 . (A. 8) 

The above relations prove that formula (A.4) really interpolates the function with appro- 
priete boundary condition 


u(nh ) = u n , n < 3 and i/(0) = H. (A.9) 

From this approximation it is straightforward to calculate the second derivative for the first 
two points: 


D 2 uq = — 9 (85^o + 66 hH — IO81/1 + 27^2 — 4^), (A.10) 

18 h z 

D 2 m = -Aj (29 u 0 + 6 hH - 54m - 2 u 3 + 27 u 2 ). (A.ll) 

18 h z 

B Supplementary material 

We have prepared a number of short movies, labelled M01 ... M13, to illustrate aspects of 
our findings, which are listed in this appendix. The movies themselves can be found at [17]. 


The first six movies show the processes depicted in figure 4 a - f: 

M01_BndryScattering_Hminus040_v020 .mov : H = —0.4, v = 0.20 < v cr (H ) 

Almost-perfect reflection of the incident antikink, which has insufficient energy to get 
over the saddle-point potential barrier. 

M02_BndryScattering_Hminus040_v0333 .mov : H — —0.4, v = 0.333 = v cr (H) 

Antikink incident at the critical velocity, leading to the creation of the saddle-point 
configuration. Note, this movie (and the associated figure 4b) is somewhat idealised, 
as in practice it is impossible to tune the initial velocity finely enough to hit the true 
critical velocity precisely. Instead, we patched together an animation up to t = 40 
with the static solution thereafter. 
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M03_BndryScattering_Hminus040_v040 .mov : H — —0.4, v = 0.40 > v cr (H) 

Single bounce, with subsequent escape of the antikink. Note that the acceleration of 
the antikink after it surpasses the potential barrier is clearly visible. 

M04_BndryScattering_Hplus090_v035.mov : H = 0.9, v = 0.35 

Single bounce, with excitation of both the H > 0 boundary mode and the internal 
mode of the antikink. 

M05_BndryScattering_Hplus090_v037.mov : H — 0.9, v = 0.37 

Single bounce exciting the boundary mode strongly enough to induce decay of the 
metastable boundary state, creating an additional kink in the bulk. 

M06_BndryScattering_Hplus090_v039.mov : H = 0.9, v = 0.39 

A similar process to M5, but here the relative velocities of the emitted kink and antikink 
are such that a bulk oscillon is formed instead of a separated kink-antikink pair. 

The next six movies scan through a range of velocities at constant H. Movies M07 - M10 show 
‘traditional’ window formation as in the full-line case (equivalent to M09). The sharpened 
edges of the windows for H < 0, caused by the presence there of a potential barrier, are 
visible on careful comparision of M7 and M8 (for H < 0) with M9 and M10 (for H > 0). For 
movies Mil and M12, H is in the region where collision-induced boundary decay is possible. 

M07_VelocityScan_Hminus040.mov : H = —0.4 

M08_VelocityScan_Hminus020.mov : H = —0.2 

M09_VelocityScan_H000.mov : H — 0 

M10_VelocityScan_Hplus020.mov : H — 0.2 

Mll_VelocityScan_Hplus090.mov : H — 0.9 

M12_VelocityScan_Hplus095.mov : H — 0.95 

Finally, movie M13 shows the slow-then-fast relaxation of the boundary mode. Four plots 
are shown: on the left, the field values next to, and slightly further from, the boundary; on 
the right, the amplitude of the field at the boundary, and an estimate cj(0,£) := 2n/T of 
its instantaneous frequency, where T is the time between successive minima of 0(0, t). 

M13_Relaxation_Hplus08393_A030.mov : H = 0.8393, A 0 = 0.3 
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